O-A053  546  MOOffC  SCHOOL  Of  ELECTRICAL  ENOINEERINO  PHILADELPHIA  PA  f/6  17/2.1 

MOBILE  communication  NETWORKS  WITH  RANGE  RESTRICTED  CHANNELS. (U) 

DEC  75  DA  MALLEN 

UNCLASSIFIED  NL 


l'»'l 

^8m548 

i 

1 

1 

1 

1 

1 

1 

1 

■ 

■ 

■ 

1 

1 

■ 

1 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

■ 

1 

1 

1 

1 

1 

1 

END 

DATE 

niMfO 

6 -78 

DOC 

CD 


University  of  Pennsylvania 
TMi;  MOORE  ECIa'OOL  OP  ELSCTRICAL  EriGINEERIMG 


merits  for  t.he  degree  of  Master  of  Science  jn  Engineering 
(for  graduate  v.’ork  in  Computer  and  Information  Sciences). 


Philadelphia,  Pennsylvania 

^fpjDec4mm,m?s  I 


niCTlirbUTION  STATEMENT  A 

Appiovod  foi  public  release; 
Distribution  Unlimited 


D D C 


m Odd 


University  of  Pennsylvania 
THE  MOORE  SCHOOL  OF  ELECTRICAL  EHGIHEERIHG 


Title  of  thesis:  Mobile  Communic'^tion  Metworks 
f '.Vith  Range  Restricted  Channels 

AbE^.tract: 

A graph  model  of  a communication  network  is  considered 
whose  nodes  represent  the  communication  stations  and  whose 
links  represent  the  communication  channels.  The  links  are 
undirected  and  there  are  no  parallel  links  or  self  loops. 
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Chapter  1 

Related  Results  from  the  Literature 

Several  authors  have  developed  procedures  for  synthesiz- 
ing highly  reliable  or  invulnerable  networks.  The  work  of 
sonie  of  these  authors  is  discussed  briefly  below  in  order  to 
gain  an  understanding  of  the  background  of  this  paper. 

Amin  and  Hakimi  (Al-ll)  prove  that  for  a graph,  G,  of  n 

nodes  and  e links,  the  network  of  G,  denoted  rJ(G)  , is  least 

vulnerable  if  the  connectivity,  w,  is  maximum  over  all  graphs 

with  n nodes  and  e links.  They  prove  further  that  N(G)  is 

maximally  survivable  if  the  independence  number,  2(G),  is 

\ 

minimiim,  v/here  the  independence  number  is  the  cardinality  of 
the  maximum  independent  set,  and  a set  ScM,  the  set  of  nodes, 
is  an  independent  set  of  G if  no  two  nodes  in  S are  adjacent 
in  G.  Amin  and  Hakimi  then  consider  the  construction  of  G 
with  a given  connectivity  and  a minimum  independence  number 
in  order  to  increase  the  network  survivability. 

Boesch  and  Thomas  (B02  and  B03)  consider  the  connecti- 
vity and  the  cohesion,  k,  and  they  prove  that  an  optimal 
damage  resistant  network  is  one  in  v/hich  w=X=2^.  Since  any 
optimal  damage  resistant  netv/ork  must  be  regular,  they  con- 
struct such  networks  for  given  values  of  n and  e by  the  union 
of  edge  disjoint  Hamiltonian  cycles.  Boesch  and  Felzer  (BOl) 
expand  this  result  to  enumerate  classes  of  optimal  damage 
resistant  networks. 

Algorithms  for  the  construction  of  optimal  daimage  re- 
sistant networks  are  given  by  Chou  and  Frank  (CKl  and  FK5) 
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such  that  the  number  of  branch  disjoint  paths  betv/een  n^  and 
n.  is  at  least  r. where  r- • is  an  element  of  the  branch  re- 

J J*  J J 

dundancy  matrix,  R={r.  •}.  Using  the  result  of  Boesch  and 

X J 

Thomas,  Chou  and  Frank  form  Hamiltonian  cycles  of  first  order 
in  their  algorithms.  They  show  that  if  all  the  node  degrees, 
d^  for  i=l,...,n,  are  even  then  the  Hamiltonian  cycles  need 
not  be  restricted  to  first  order. 

Frank  (FKl)  generates  a class  of  networks  by  a random 

process  where  all  of  the  networks  are  symmetric  in  probability, 

develops  a model  for  an  attack  strategy,  and  considers  the 

$ 

survivability  of  the  networks  generated.  Here,  he  considers 
the  effects  of  uniform  and  'nonuniform  links,  and  links  having 
a distance  bias.  In  a later  work,  he  (FK2)  considers  a class 
of  networks  having  perfectly  reliable  links  and  nodes  having 
a probability  p of  existing.  The  survivability  of  the  network 
is  examined  for  p close  to  unity. 

In  two  more  recent  articles,  Frank  (FK3  and  FK4)  examines 
command  and  control  communication  systems.  His  network  analy- 
sis begins  v.ith  modeling  the  physical  environment  and  ends 
with  the  solution  of  the  analysis  problems.  Several  measures 
of  survivability  are  considered,  including  the  number  of  nodes 
in  the  largest  communicating  section  of  the  network  after  the 
attack. 

Fratta,  et  al,  (FRl)  consider  the  design  of  networks  in 
which  the  connection  probability  for  the  terminal  node  pairs 
in  the  worst  conditions  exceed  some  predetermined  value.  In 
their  study  the  nodes  of  the  network  are  complc  tely  reliable 
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while  the  links 


are  nonuniforrnly  unreliable. 


Hansler  (ILll)  x^^'^sents  an  algorithm  that  recursively 
calculates  the  probability  that  all  paths  between  two  nodes 
in  a communication  netv/ork  have  been  interrupted.  His  algo- 
rithm applies  to  networks  v/ith  nonuniforrnly  unreliable  links 


(HAk) 


only  unreliable  links  and  nodes,  but  networks  whose  reliabil- 
ity also  varies  with  the  netv/ork  topology.  Their  approach  is 
to  increase  the  reliability  of  the  netv.-ork  by  removing  unre- 
liable links  and  adding  more  reliable  ones  of  lower  cost. 

This  approach  follows  that  of  Steiglitz,  et  al,  (STl)  who 
also  develop  a procedure  for  designing  optimally  reliable 
networks  subject  to  a cost  constraint. 

Vlilkov  gives  an  approximation  method  for  generating 

nonregular  degree,  maximally  connected  networks  of  minimum 
diameter,  where  the  diameter  is  the  maximum  length  of  any 
shortest  path  in  the  graph,  l.'illiams  (b'Ll)  cites  three  ways 
to  enhance  the  survivability  of  a network;  hardening  the  ele- 
ments, adjusting  the  location  of  the  elements,  and  adding  re- 
dundant elements.  In  his  consideration  the  choice  of  methods 


is  dictated  by  the  costs  involved,  and  he  includes  a procedure 
for  approximating  such  costs. 

The  presentation  above  is  by  no  means  an  exhaustive  sum- 
mary of  the  literature  dealing  with  survivable  netv/orks.  The 
reader  is  directed  to  Wilkov  (oKl)  for  an  excellent  summation 
of  a wide  range  of  authors'  works,  h'ilkov  also  includes  an 
oxten.sive  bibliography. 


Chapter  2 
Mobile  Networks 

2.1,:  Introduction 

In  this  paper  a class  of  networks  termed  "mobile  net- 
v/orks"  is  studied.  A mobile  network  is  defined  as  a network 
whose  nodes  have  positions  which  vary  with  time.  The  mobile 
networks  studied  are  characterized  by  the  following  assump- 
tions! the  links  of  the  network  are  undirected;  the  network 
has  no  parallel  links;  the  network  has  no  self  loops;  each 
node  of  the  netw'ork  has  a given  initial  and  terminal  point. 
2.2i  Problem  statement 

The  problem  dealt  with  here  can  be  divided  into  two 
major  parts.  First,  develop  a procedure  which  moves  the 
network  from  its  given  initial  position  to  its  given  terminal 
position  subject  to  the  following  restrictions! 

1)  The  links  have  a fixed  maximum  length. 

2)  The  nodes  may  move  no  more  than  a fixed  maximum 
distance  during  one  iteration. 

3)  Only  a fixed  maximum  percentage  of  the  nodes 
may  move  at  one  time. 

4)  The  nodes  may  not  assume  intermediate  positions 
coincident  with  other  nodes. 

The  second  part  of  the  problem  is  to  determine  the  sur- 
vivability of  the  network  at  its  initial  and  terminal  position 
and  at  each  intermediate  position.  For  this  analysis  both 
the  links  and  the  nodes  are  assumed  to  be  nonunifonnly  unre- 
liable. The  reliability  of  a node  is  inversely  proportional 


4 


to  its  degree,  and  the  reliability  of  a link  is  inversely 
proportional  to  its  length.  The  network  survivability  v/ill 
be  determined  by  subjecting  the  network  to  a series  of  random 
attacks  and  computing  the  percentage  of  these  attacks  in  which 
all  of  the  surviving  nodes  remain  in  one  connected  component. 
2.3;  Overview 

This  paper  presents  a solution  to  the  problem  stated 
above.  Chapter  3 describes  the  algorithm  this  author  devel- 
oped v/hich  successfully  translates  the  mobile  network  from 
the  initial  .to  the  terrunal  position  while  calculating  the 
network  survivability 'after  each  intermediate  movement. 

Chapter  4 gives  a detailed  explanation  of  the  portions  of  the 
algorithm  which  perform  the  actual  movement  of  the  nodes. 

The  criteria  used  to  deterTnine  which  i.odes  to  move  is  also 
explained.  Chapter  5 describes  the  representative  networks 
chosen  to  test  the  algorithm.  For  each  such  netv/ork,  upper 
and  lower  bounds  on  the  link  and  node  reliability  are  calcu- 
lated and  shown  in  Chapter  5*  Chapter  6 gives  the  results 
of  the  survivability  an-xlysis  for  the  networks  tested  and 
the  conclusions  drawn  from  this  study. 

2.4t  An  application  jf  mobile  netv/orks 

Vihen  using  a computer  as  an  aid  in  solving  such  a netv/ork 
problem,  it  is  understood  that  the  original  problem  has  been 
reduced  to  a graph  model  of  the  netv/ork.  Indeed,  the  problem 
statement  of  Section  2.2  mentions  only  "links"  sind  "nodes" 
and  the  attributes  they  are  given  in  this  particular  problem. 
Thus  the  author  has  presented  a graph  model--but  of  v/hat? 
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VJhile  there  are  many  possible  applications  of  mobile  netv.'orks, 
only  one  will  be  presented  here. 

Assume  that  a number  of  military  units  are  advancing 
toward  their  enemy.  The  enemy  is  believed  to  have  no  prior 
knowledge  of  the  unit  locations  although  he  undoubtably  has 
a limited  ability  to  detect  these  positions.  During  the  ad- 
vance, it  is  vitally  important  that  the  units  be  able  to 
communicate  with  each  other,  either  directly  or  by  relaying 
messages . 

This  situation  is  perfectly  suited  for  the  mobile  network 
model.  Each  unit,  or  station,  in  a communication  network  can 
be  represented  as  a node  in  the  graph.  The  communication 
channel  joining  two  stations  in  the  netv.ork  can  be  represented 
by  a link  in  the  graph. 

One  restriction  on  the  mobile  netv/orks  imposed  by  the 
problem  statement  is  that  the  links  have  a fixed  maximum 
length.  In  Chapter  1 reference  was  made  to  Frank  (FKl)  who 
discusses  networks  whose  links  have  a distance  bias.  Both 
of  these  concepts  readily  fit  the  military  communication  net- 
v/ork  application.  The  maximum  length  corresponds  to  the  max- 
imum transmission  distance  of  the  transmitter  at  each  station. 
The  distance  bias  corresponds  to  the  lessening  of  reliability 
of  the  transmitted  signal  as  its  distance  from  the  source  in- 
creases. This  lessening  of  reliability  results  from  signal 
v/eakness,  terrain  interference,  enemy  jamming,  etc.  Relating 
these  tv/o  concepts,  the  maximum  transmission  distance  is  de- 
fined as  the  range  corresponding  to  the  signal  threshold 
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reliability.  Beyond  this  range'  (or  below  this  threshold)  the 
coiiimuni cation  quality  via  the  channel  is  considered  to  be 
unacceptable . 

The  nodes  of  the  network,  ratho than  being  perfectly 
reliable,  have  a survivability  which  is  inversely  proportional 
to  their  degree.  P’or  the  communication  network,  the  degree 
of  a node  represents  the  number  of  units  v/ithin  the  maximum 
transmission  range  of  a given  unit.  As  the  number  of  units 
v/ithin  this  range  increases,  the  enemy's  ability  to  detect 
the  unit's  location  also  increases — due  not  only  to  the  in- 
crease in  his  target  density  but  also  to  the  probable  increase 

in  message  traffic.  Therefore,  as  the  enemy's  probability  of 

\ 

success  increases,  the  survivability  of  the  unit  decreases. 
Because  of  this  increase  in  vulnerability  of  densely  populated 
areas,  the  restriction  that  no  tv/o  nodes  in  the  network  model 
may  be  coincident  applies  to  the  communication  network.  By 
extending  the  restriction,  a minimum  separation  distance  is 
established  v/hich  helps  position  the  nodes  over  a larger  area, 
thus  decreasing  the  density  of  the  corresponding  military 
units  and  increasing  their  survivability. 

The  netv/ork  survivability  criterion  used  here  also  ap- 
plies to  the  military  communication  netv/ork.  After  an  attack, 
the  surviving  units  must  be  able  to  react  to  orders  initiated 
by  the  commander.  Autonomous  units — those  units  v/hich  sur- 
vive but  cannot  communicate  with  the  commander — reduce  his 
effectiveness  since  he  is  forced  to  assume  that  those  units 
v/ore  also  destroyed. 
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TKe  N«  I' 


3.1:  The  ^ ^ j j . 

The  alcori  if  * , r 

algorithm  from  tr.«  . , t.-  ' 
computed  by  deti  rr  i -,; 
v/ork  after  it  bar.  l*-*--  ^ 

procedure  ia  accumt* . * -r.. 

(FK4).  Frank  devt  . : « 
the  number  of  cor.p«  •.* 
and  nodes  are  cut  , « f • ’• 
the  target  of  a rand;  . 
mobile  networks  cltar. 

One  modi  flee  ’ ’ r ■ 

After  determini  »<r.c^  " i- 
tions  are  calculated  tor 
temporary  position  :• 

nected.  This  check  ;e  *.  • r », 
Algorithm**  and  bypast.  1:.^ 

The  bypass  is  acconpl  lar.f  c Z'- r .4 
which,  when  Off  ( FAUi:.  < , <i.r»c*. 
the  number  of  component:  * i*  - 

On  (TRUE),  directi:  th»  pr.<»a-' 
components  in  a protati..  . 
random  attack. 

3.2:  S^/mbolorv 

A listing  of  ‘.he  ; . 


Algorithm  appears  in  Section  3*3*  A flowchart  of  the  proce- 
dure is  given  in  Figure  I.  The  following  symbol  definitions 
apply  to  symbols  found  both  in  the  algorithm  and  in  the 
flowchart. 

n the  number  of  nodes  in  the  network 

the  remaining  distance  each  node  must  travel 
Fj^  terminal  netw'ork  position 

dp  minimum  node  degree  in  Fj^ 

initial  network  position 
dj  minimum  node  degree  in 

Tj^  temporary  netv/ork  position 
d^  minimum  node  degree  in  Tj^ 
dj^  minimum  net\vork  degree 

Cj^  current  network  position 

Rq  reliability  of  Cj^ 

r number  of  random  attacks 

B number  of  random  attacks  survives 

nt  maximum  number  of  nodes  moving  per  iteration 
nodes  that  move  in  current  iteration 
MA^  unmoved  nodes  not  in  main  component 
MBg  moved  nodes  not  in  main  component 
Vnv  nodes  not  in  main  component  after  move 
t time 

t maximum  allowable  time 

max 

hfj^  multiplication  factor 

mf  maximum  value  of  hf 

D.  . distance  between  nodes  i and  j 

i J 
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rn  maximum  transmission  range 

md  maximum  movement  distance 

a minimum  separation  distance 

“B.lt  The  network  translation  algorithm 

0.  Preliminary  calculations. 

a)  Calculate  for  i=l,...,n. 

b)  If  Fj^  is  not  connected,  HALT.  Calculate  dp. 

c)  If  Ij^  is  not  connected,  HALT.  Calculate  dj . 

d)  dj^  = mint  dp,  dj}. 

e)  Set  t = 0 and  Cj^  = Ij^. 

1 . Random  attack . 

a)  Set  s = 0.  Conduct  r random  attacks  on  Cj^.  If  after 
the  ith  attack  the  number  of  components  in  Cj^  is  one, 
set  s = s + 1. 

b)  “ r*  " 

2.  Node  movement. 

a)  If  Cj^  = Fj^,  HALT. 

b)  hf^  = 0 for  i=l,...,n. 

c)  Let  M^.  be  the  nodes  v/ith  the  largest  D^. 

^ n ) 

d)  Calculate  T^^  by  moving  a distance  2 md 

closer  to  Fj^  for  i=l,...,nt. 

e)  Insure  that  no  nodes  in  Tj^  are  closer  than  a. 

3.  Temporary  position  checkout. 

a)  Is  Tj^  connected?  If  so,  set  ic  = 0 and  go  to  3b). 

1 ) Separate  V into  and  MBq . 

nv  a p 

2)  Set  RO  . = M . - V . Sort  d.  for  i € MO  in 

nt-nv  nt  nv  i 
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descending  order.  Let  no  = nt  - nv. 

3)  Do  iteratively  for  i=l,...,nv.  If  i C MB  go  to 
b]  below, 

a]  i € MA: 

l]  j = where  k = no  - i . 

3]  \ = MA^. 

if]  hf,  =0. 

5J  Move  i 2 md  closer  to  F^.  Store  the 

result  in  T^. 

6]  Set  ic  = 0 and  go  to  3b) . 

b]  i € MB: 

» 

1]  hf,  = hf,  +1. 

_ ^ ^ (-hf.) 

2J  Move  i closer  to  by  2 md.  Store 

the  result  in  . 

3]  Set  ic  = 0. 

b)  If  d^  > dj^  go  to  4. 

1)  ic  = ic  + 1 . 

2)  Let  j = nt  - ic,  and  set  T.  = C.. 

J J 

3)  Let  i be  such  that  d^  < dj^.  If  i € go  to  b], 

a]  Move  i closer  to  F^^  by  2^”^^i^md.  Store  the 

result  in  T^,  and  go  to  3b). 

b]  hf^  = hf^  +1.  If  hfj^  < mh,  go  to  a]  above. 

1]  Move  i from  closer  to  F^  by  l-iiiJ.  Store 

the  result  in  . 

2]  Locate  a node  k such  that  if  < rn 
then  d,  > d.,. 

1 ■“  iN 
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3]  Move  i toward  k until  D.,  < rn. 

XK 


4]  Go  to  3b) . 

4.  Preparation  for  the  next  time  step. 

a)  t = t + 1.  i 

b)  Set  Cj^  = 

c)  Go  to  1 . 


12 


f 


b'etY/ork  Translation  Algorithm  Flowchart 


START , 


t = 0 


Calculate 
i=l n 


-Are 

.connected?^ 


{HAIT 


Calculate  d^, 

dj^  = min[dp,dj] 


Conduct  random 
attacks  on  Cj^. . 

Calculate  R« 


Program  halts: 

0 : normal 

termi nati on 

1 : end  position 
not  connected 

2:  maximum  run 
time  exceeded 


C = F 

N- 


Mo 
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Chapter  4 
Node  Movement 

4.1 t Inbroduction 

The  most  important  steps  of  the  net\';ork  translation  al- 
gorithm— and  the  most  difficult  steps  to  comprehend--are  the 
ones  dealing  v/ith  the  movement  of  nodes  and  the  checkout  of 
the  temporary  positions  generated.  Reference  was  made  in 
Ch  pter  1 to  Williams ' (iVtl)  three  methods  for  enhancing 
network  survivability.  Recall  that  the  methods  are: 

1)  Physically  harden  the  network  elements. 

2)  Adjust  the  location  of  the  network  elements. 

3)  Add  redundant  network  elements. 

The  addition  of  redundant  network  elements  is  considered  an 
impractical  means  of  improving  network  survivability  for  the 
mobile  netv.'orks  studied  here.  This  chapter  will  describe 
how  the  algorithm  adjusts  the  location  of  the  network  ele- 
ments to  improve  survivability,  and  Chapter  5 v/ill  deal  with 
the  problem  of  hardening  the  network  elements. 

Amin  and  Hakimi  (Ai'Il)  proved  that  the  survivability  of 
a network  can  be  increased  by  increasing  the  connectivity. 

The  connectivity  is  increased  by  adding  links  to  those  nodes 
of  minimum  degree.  Unfortunately,  in  mobile  networks  with 
range  restricted  links,  one  cannot  simply  add  links  to  in- 
crease the  connectivity  since  all  of  the  links  meeting  the 
range  restriction  already  appear  in  the  netv/ork.  This  author 
did  use  the  approach  of  Amin  and  Hakimi,  however,  in  conjunc- 
tion v/ith  Vi/illiaras'  proposal  to  adjust  the  location  of  the 
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network  elements.  As  a result,  the  survivability  of  the  mo- 
bile netv/orks  tested  is  increased  by  moving  the  nodes  in  such 
a way  as  to  increase  the  netv/ork  connectivity.  Since  the 
mobile  networks  studied  have  fixed  initial  and  terminal  posi- 
tions, the  connectivity  of  each  of  these  end  positions  can 
be  easily  determined.  These  two  connectivities  are  then  com- 
pared, and  their  minimum  value  is  defined  to  be  the  minimum 
netv/ork  connectivity,  or  simply  the  network  connectivity. 

Thus,  at  the  intermediate  positions — those  positions  where  the 
location  of  the  elements  can  be  varied--each  node  is  required 
to  have  a degree  at  least  as  great  as  the  network  connectivity. 
4.2 » Choosing  v/hich  nodes  to  move 

I 

The  network  translation  algorithm  was  formulated  in  such 
a v/ay  that  it  attempts  to  move  the  netv/ork  from  its  initial 
position  to  its  terminal  position  in  the  minimum  number  of 
iterations.  In  this  context,  an  iteration  is  a movement  of 
nt  nodes  where  nt  corresponds  to  the  maximum  percentage  of 
nodes  allowed  to  move  concurrently.  Since  the  most  direct 
route  from  one  point  to  another  is  simply  a straight  line, 
the  algorithm  will,  when  possible,  adhere  to  straight  line 
motion  for  each  node  moved.  As  v/ill  be  discussed  later,  there 
are  two  possible  situations  v/hich  result  in  nonlinear  motion. 

The  basic  criterion  for  selecting  which  nodes  to  move  is 
the  distance  between  a node's  current  position  and  its  terminal 
position.  This  distance  is  calculated  by 


v/hcre  d is  the  number  of  dimensions  and  p and  q are  the  points 
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betv/een  v/hich  the  distance  is  being  measured.  The  algorithm 
locates  those  nt  nodes  v;ith  the  greatest  remaining  distance, 

D , and  translates  them  the  maximum  allowable  movement  dis- 
tance , md,  along  the  line  connecting  their  current  position 
and  their  terminal  position. 

The  calculated  positions  of  the  nodes  moved,  along  with 

I 

the  current ■ positions  of  the  nodes  not  moved,  represent  the 
temporary  netv/ork  position.  Tv/o  tests  are  conducted  on  the 
network's  temporary  position:  (a)  is  it  in  one  connected 
component;  and  (b)  is  the  degree  of  each  node  at  least  as 
great  as  the  network  connectivity.  Sections  4.3  and  4.4  ex- 
plain the  procedures  used  when  these  tests  fail. 

4.3:  GeneratinfT:  a connected  comnonent 

In  a mobile  net^vork  v/here  less  than  100^  of  the  nodes 
move  in  each  iteration,  it  is  quite  possible  that  some  motion 
will  result  in  disconnecting  the  network.  There  are  three 
cases  that  v/ill  dissolve  the  network  into  more  than  one  com- 
ponent: 

Case  1:  A portion  of  the  moving  nodes  move  beyond 
the  range  of  the  main  component. 

Case  2:  A portion  of  the  nonmoving  nodes  is 

"stranded"  as  the  main  component  moves 
away  from  it. 

Case  3:  Some  combination  of  Case  1 and  Case  2, 

The  discussion  below  will  focus  first  on  a Case  1 node  and 
then  on  a Case  2 node.  Discussing  Case  3 would  be  redundant 
since  it  employs  the  techniques  of  both  Case  1 and  Case  2. 


The  distance  a node  moves  during  an  iteration  is  given 
by  2'~  Vd , v/here  md  is  the  maximum  allowable  movement  dis- 
tance and  hf  is  a multiplication  factor.  Initially,  hf  is 
set  to  zero  for  each  node  and  the  moving  nodes  are  translated 
a distance  2^md,  or  simply  md.  After  the  initial  movement, 
if  the  netv/ork  is  no  longer  connected,  hf  is  incremented  by 
one  for  each  Case  1 node  and  the  temporary  position  of  each 
such  node  is  recomputed.  Since  the  Case  1 node  has  moved  be- 
yond the  main  component,  halving  the  distance  it  moves  tends 
to  bring  it  back  into  the  range  of  the  main  component.  The 
halving  process  continues  a finite  number  of  times — a number 

determined  by  the  user  based  on  the  maximum  movement  distance 

\ 

and  the  maximum  range  of  the  links. 

A Case  2 node  presents  an  entirely  different  problem. 

The  nodes  that  had  been  adjacent  to  the  Case  2 node  connecting 
it  with  the  main  component  have  moved  beyond  its  link  range. 
The  problem  would  be  simply  resolved  by  moving  the  Case  2 
node;  however  moving  it  will  increase  the  number  of  nodes 
moved  to  a value  greater  than  the  maximum  allowed.  Therefore 
some  node  v/hich  did  move  must  be  moved  back  to  the  position 
it  maintained  at  the  start  of  the  iteration  before  the  Case 
2 node  can  move.  A criterion  for  selecting  the  node  that 
must  move  back  is  needed.  Consider  the  degree  of  the  moved 
nodes  I a node  of  high  degree  has  a lov/er  likelihood  of  being 
in  a "critical"  position — that  is,  a position  v/hich  makes  it 
the  lone  connector  between  some  other  node  and  the  main  com- 
ponent. Using  this  reasoning,  the  nctworl:  translation 
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algorithm  finds  the  node  of  highest  degree  among  the  moved 
nodes  and  returns  it  to  the  position  it  held  at  the  beginning 
of  the  iteration.  The  Case  2 node  is  then  moved  forward  and 
the  maximum  percentage  of  nodes  allowed  to  move  is  not  ex- 
ceeded. 

Each  time  there  is  more  than  one  component  in  the  net- 
work, the  Case  1 nodes  are  repositioned  to  place  them  in  the 
main  component  before  the  Case  2 nodes  are  considered.  One 
node  is  brought  back  into,  or  moved  forward  into,  the  main 
component  at  a time,  and  after  each  such  adjustment  the  net- 
work is  reevaluated  i-n  order  to  prevent  unnecessary  changes 

since  correcting  one  node  may  have  the  side  effect  of  cor- 

> 

recting  others. 

U.U:  Increasin/r  the  connectivity 

After  the  algorithm  verifies  that  the  network  is  com- 

fi 

posed  of  one  connected  component,  the  problem  of  connectivity 
of  individual  nodes  is  studied.  If  there  is  a node  whose  de- 
gree is  less  than  the  network  connectivity,  the  position  of 
this  node  v/ill  be  altered  to  increase  its  degree. 

The  node  whose  connectivity  must  be  increased  will  move 
from  its  temporary  position  back  to\/ard  the  position  it  main- 
tained at  the  start  of  the  iteration,  moving  ^ unit  each  step. 
After  each  such  move,  the  degree  of  the  node  v/ill  be  recalcu- 
lated and  the  procedure  v/ill  end  if  the  degree  is  at  least  as 
great  as  the  network  connectivity.  If  the  degree  requirement 
is  still  not  satisfied,  the  node  continues  moving  backward  in 
i unit  steps.  The  process  stops  if  the  node  reaches  a point 
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only  J unit  from  where  it  began  the  iteration.  Having  ex- 
hausted the  possibility  of  each  half  unit  interval  v/ithin 
the  maximum  movement  range  along  the  line  connecting  its  start 
point  and  its  terminal  point,  the  node  enters  a routine  which 
results  in  nonlinear  motion  in  order  to  satisfy  the  degree 
requirement  on  the  node. 

The  node  is  moved  from  its  starting  point  for  the  itera- 
tion a distance  -Jrnd  toward  its  terminal  position.  From  this 
position  it  searches  for  the  closest  node  that  would  raise 
its  degree  to  at  least  the  minimum  netv/ork  connectivity.  When 
such  a node  is  located,  the  moving  node  is  translated  along 
the  line  connecting  its  current  position  and  the  node  found 
until  the  two  nodes  are  separated  by  a distance  less  than  or 
equal  to  the  maximum,  link  range. 

The  discussion  above  assumes  that  the  node  v/hose  degree 

<r 

must  be  increased  is  one  of  the  nodes  moved  during  the  itera- 
tion. If  in  fact  it  was  not  moved,  it  is  treated  as  a.  Case  2 
node,  in  the  terminology  defined  in  Section  4.3,  and  some 
other  node  is  returned  to  its  starting  point  so  the  Case  2 
node  can  be  moved  to  increase  its  degree. 

4.5t  I^iSintaining  a minimum  senaration  distance 

Each  time  a node  is  moved,  its  distance  from  its  neigh- 
bors is  checked.  If  the  distance  between  the  moved  node  and 
any  other  node  is  less  than  the  minimum  allowable  separation, 
the  position  of  the  moved  node  must  be  varied.  The  variation 
in  movement  resulting  from  this  situation  is  the  second  in- 
stance of  nonlinear  motion  contained  within  the  network 
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translation  al^ox 


The  moved  node  » : i 

neichbor  alone  the  li  ■ 
the  distance  betv.et-fi  'f 
separation  distar.ee. 
in  position  will 
node  exceedinf^  the  aia 
distance  moved  is  recr • • 
has  been  exceeded,  th* 
point  until  both  the  • . 

minimum  separatior  re: 


Chapter  5 

Mobile  Network  Reliability 

5.1;  Hardening:  net'v/ork  elements 

One  method  cited  by  b'illiarns  (WLl)  for  increasing  the 

network  survivability  is  hardening  the  network  elements. 

There  are  two  basic  ways  to  simulate  the  hardening  of  netv/ork 

elements  for  the  purpose  of  network  analysis.  The  more  subtle 

approach  is  given  by  Frank  (FKl).  In  his  netv/ork,  station  s 

is  destroyed  if  it  is  hit  by  K v.'eapons.  Similarly,  channel  c 

s 

is  destroyed  if  it  is  hit  by  weapons.  Frank  assumes  that 
each  weapon  can  produce  at  most  one  hit,  and  the  probability 

I 

of  a network  element  being  hit  by  any  given  v;eapon  is  inde- 
pendent of  the  number  of  hits  the  element  has  already  sus- 
tained. This  approach  relates  the  hardness  of  an  element  to 
the  number  of  v^eapon  hits  required  to  totally  destroy  it. 

The  effect  of  varying  degrees  of  hardness  can  be  compared  by 

changing  the  values  of  K and  K . Such  an  approach  would  be 

s c 

quite  realistic  for  a military  application  since  in  a conven- 
tional warfare  situation  a unit  would  probably  not  be  totally 
destroyed  by  a single  weapon  hit. 

The  second  way  to  simulate  hardening  the  netv/ork  elements 
is  the  Monte  Carlo  method  used  by  Baran  (BAl).  Baran  tests 
the  survivability  of  networks  after  establishing  values  for 
the  node  and  linlc  reliabilities.  By  varying  these  reliabil- 
ities various  degrees  of  element  hardness  can  be  compared. 

Despite  the  appropriateness  of  Frank's  method  to  the  ap- 
plication developed  in  Chapter  2,  Baran 's  approach  is  used  ir 
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this  paper.  Franlc's  method  is  viewed  as  being  too  cumbersome 
for  the  analysis  included  here. 

5.2t  Node  reliability 

In  the  netv;ork  translation  algorithm,  the  survival  pro- 
bability of  node  i is  given  by 

.(100.0  - t„)io"^ 

= 1.0  - [ n - 1 ^^i 

v;here  t^  is  the  node  threshold  reliability  and  d^  is  the  de- 
gree of  node  i. 

Upper  and  lower  bounds  on  the  node  reliability  can  be 
established.  Since  the  reliability  of  the  node  is  inversely 
proportional  to  its  degree,  the  lower  bound  on  the  node  relia- 
bility corresponds  to  a node  of  maximum  degree.  In  the  mobile 
networks  considered  here,  where  there  are  no  parallel  links 
and  no  self  loops,  the  maximum  degree  of  any  node  in  a net- 
v/ork  of  n nodes  is  n-1 . The  given  value  of  the  node  thresh- 
old reliability  is  made  to  correspond  to  a degree  of  n-1.  An 
upper  bound  on  the  node  reliability  corresponds  to  a node  of 
minimum  degree.  As  discussed  earlier,  the  minimum  degree  of 
any  node  is  the  calculated  value  of  the  network  connectivity. 

The  algorithm  v/as  tested  on  the  tv/o  networks  that  will  be 
presented  in  Chapter  6.  The  smaller  network  has  8 nodes  and 
the  larger  netv/ork  has  25  nodes.  In  each  network,  the  minimum 
network  connectivity  is  three.  Table  I gives  upper  limit 
values  for  the  node  survival  probability  for  each  netv/ork 
where  the  lov/er  limit  is  the  threshold  reliability. 

5.3i  Link  reliability 

In  the  netv/crk  translation  algorithm,  the  survival  pro- 
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Bounds  on  the  Node  Reliability 
Lower  Upper 


n = 8 

n = 25 

70.0% 

87.14^5 

96.25% 

80.0 

91.43 

97.50 

90.0 

95.71 

98.75 

100.0 

100.00 

100.00 

Table  I 


bability  of  link  j is  given  by 


>-2 


(100.0  - t,)l0 

Pj  = - C ^“ilc 

v/here  is  the  link  threshold  reliability,'  rn  is  the  maximum 
link  range,  and  is  the  distance  between  nodes  i and  k, 
the  endpoints  of  the  jth  link. 

Upper  and  lower  bounds  on  the  link  reliability  can  also 
be  established.  Since  the  reliability  of  the  link  is  inver- 
sely proportional  to  its  length,  the  lower  bound  on  link  re- 
liability corresponds  to  a link  of  maximum  length.  The  upper 
bound  on  the  link  reliability  corresponds  to  a link  of  minimum 
length  v/here  the  minimum  length  is  given  by  the  minimum  sepa- 
ration distance.  For  the  tv/o  networks  tested  the  maximum 
range  was  the  same,  as  was  the  minimum  separation.  Table  II 
gives  upper  limit  values  for  the  link  survival  probability 
for  varying  lo\.’er  limits,  where  the  lov/er  limit  is  equal  to 
the  given  value  of  the  threshold  reliability. 


ii 


25 


Bounds  on  the  Link  Reliability 


Lower 

Upoer 

50.0% 

j 

93.75%' 

60.0 

95.00 

70.0 

96.25 

80.0 

97.50 

90.0 

98.75 

100.0 

100.00 

Table  II 
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Chapter  6 

Results  and  Conclusions 

6.1;  The  mobile  networks  tested 

Two  mobile  netv/orks  v/ere  selected  at  random  to  check  the 
network  translation  algorithm.  Figure  II  shows  the  graph 
model  of  the  initial  and  terminal  positions  of  the  mobile 
network  of  8 nodes.  The  algorithm  translated  the  netv/ork 
from  its  initial  to  its  terminal  position  in  five  iterations 
v/hich  is  the  minimum  number  possible  for  the  given  parameters. 

Figure  III  shows  the  graph  model  of  the  initial  and  ter- 
minal positions  for  the  mobile  network  of  25  nodes.  The  algo- 
rithm completed  the  translation  of  the  network  in  13  iterations. 
The  minimum  number,  of  iterations  for  the  netv/ork  under  the 
given  parameters  is  12;  however,  only  two  nodes  remained  to 
be  translated  after  12  iterations  and  the  greater  remaining 
distance  of  these  two  nodes  v/as  O.698  units. 

6.2}  A standard  for  comparison  of  results 

The  network  survivability  criterion  used  to  measure  the 
survivability  of  mobile  netv/orks  is  similar  to  that  used  by 
Baran  (BAl)  in  his  analysis  of  distributed  communication  net- 
works. He  judged  the  survivability  of  a network  by  the  per- 
centage of  nodes  in  connection  with  the  largest  single  group 
of  surviving  nodes,  and  he  considered  small  groups  of  isolated 
nodes  to  be  ineffective.  Baran  showed,  using  Monte  Carlo 
simulation,  that  a highly  reliable  network  can  be  constructed 
using  unreliable  links  when  a satisfactory  redundancy  level 
is  chosen.  He  defines  the  redundancy  level  as  the  link-to-node 


27 


ratio  in  an  infinite  size  array.  Baran  proved  that  for  per- 
fectly reliable  nodes,  links  of  only  50^  reliability  cause 
little  system  degradation  for  redundancy  three  or  four. 

Because  of  the  similarities  between  the  mobile  networks 
tested  and  Baran' s distributed  netv/orks,  it  is  possible  to 
use  his  results  to  measure  the  accuracy  of  the  results  gener- 
ated by  the  network  translation  algorithm.  However,  before 
proceeding  an  understanding  of  the  differences  in  the  two 
approaches  must  be  gained.  From  Baran* s definition  of  the 
redundancy  level  it  is  clear  that  the  redundancy  level  of  a 
node  is  twice  the  node's  degree.  For  networks  with  fixed 
initial  and  terminal  positions  and  range  restricted'ilinks 
the  connectivity  can  be  determined  but  not  increased.  Only 
the  connectivity  at  intermediate  positions  can  be  varied,  and 

it  cannot  be  rea'sonably  maintained  at  a higher  level  than  the 

<!• 

minimum  connectivity  of  the  two  end  positions.  As  discussed 
earlier  this  procedure  is  used  to  establish  the  network  con- 
nectivity. If  the  network  connectivity  is  six  or  eight — 
corresponding  to  a redundancy  level  of  three  or  four — then 
the  result  Baran  proved  concerning  systems  with  links  of  5^% 
reliability  could  be  checked.  However,  for  a network  connecti- 
vity of  less  than  six,  a higher  degree  of  degradation  must  be 
expected.  The  netv/orks  evaluated  by  the  network  translation 
algorithm. both  have  a network  connectivity  of  three.  V/hen  the 
results  of  the  survivability  tests  on  these  networks  are  shown 
in  Section  6.3  they  will  be  compared  to  Baran' s equivalent 
network  with  a redundancy  level  of  1.5* 
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An, Eight  Node  Mobile  Ketv/orli 


Terminal  Initial 

Position  Position 


Number  of  nodes  moving  per  iteration;  4 

Maximum  link  range ; 4.0 

Maximum  movement  distance;  4,0 
Minimum  network  connectivity;  3 


Figure  II 
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A Tv/enty-five  Ilode  Mobile  Netv;ork 
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Parameter  Values: 

Number  of  nodes  moving  per  iteration:  12 

Maximum  link  range:  4.0 

Maximum  movement  distance:  4,0 

Minimum  network  connectivity:  3 


Figure  III 
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The  v/ay  in  which  the  link  and  node  reliabilities  are  de- 
termined in  the  network  translation  algorithm  will  also  cause 
the  results  generated  to  deviate  from  Baran's.  If  the  link 
threshold  reliability  is  50%,  only  those  links  of  maximum 
length  will  be  50'%  reliable.  Clearly,  not  all  of  the  links 
will  be  of  maximum  length.  The  links  of  less  than  maximum 
length  will  be  more  tha.n  50%  reliable , thus  tending  to  reduce 
the  network  degradation.  Similarly  for  any  value  of  the  node 
threshold  reliability,  only  those  nodes  of  maximum  degree 
will  exhibit  the  threshold  reliability.  Nodes  of  lower  de- 
gree will,  therefore,-  also  reduce  the  network  degradation. 

6.38  Mobile  network  reliability  results 

Figures  IV  and  V are  graphs  of  the  mobile  network  survi- 
vability generated  for  varying  values  of  link  and  node  relia- 
bility, respectively.  The  results  are  shown  for  both  the  8 

G 

and  the  25  node  networks  as  continuous  functions.  Clearly, 
these  functional  representations  are  broad  extensions  of  a 
relatively  small  number  of  discrete  samples.  The  continuous 
curves  serve  as  a visual  aid  but  have  not  been  rigorously  de- 
termined. In  each  figure,  the  curve  labeled  "DCN"  corresponds 
to  Baran's  results  for  his  Distributed  Communication  Netv^orks. 

The  mobile  network  results  shown  in  Figures  IV  and  V are 
averages  for  each  network.  For  given  values  of  the  link  and 
node  reliabilities,  each  network  at  each  iteration  underwent 
a series  of  50  random  attacks.  The  plotted  curves  correspond 
to  the  average  calculated  value  of  the  network  survivability 
over  all  iterations. 
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, Network  Survivability 
As  a Function  of  Link  Reliability 


Figure  IV 
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Network  Survivability 
As  a Function  of  Node  Reliability 
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6. hi  Concluaionr; 

This  paper  presented  an  algorithm,  the  netv/ork  transla- 
tion algorithm,  v;hich  successfully  translates  a mobile  network 
from  a fixed  initial  position  to  a fixed  terminal  position 
while  maintaining  an  established  minimum  netv/ork  connectivity 
at  each  intermediate  iteration.  Tv/o  representative  sample 
networks  were  chosen  to  test  the  algoritlim.  Based  on  the  re- 
sults generated  by  these  tests,  the  algorithm  is  felt  to  be 
very  efficient  in  so  far  as  the  procedure  progresses  to  a 
conclusion  in  an  optimal  number  of  iterations. 

The  reliability  of  the  mobile  networks  closely  parallels 
the  reliability  of  the  distributed  networks  of  equivalent  re- 
dundancy studied  by  Baran  (BAl),  as  evidenced  by  the  surviva- 
bility curves  in  Figures  IV  and  V.  This  similarity  was  ex- 
pected since  the  requirement  of  a minimum  degree  on  each  node 
is  the  foundation  of  the  definition  of  distributed  netv/orks. 

It  is  also  evident  from  Figure  V that,  for  a small  mobile 
netv/ork,  a decrease  in  the  node  survivability  docs  not  affect 
the  network  survivability  as  much  as  does  a cori'esponding  re- 
duction of  the  node  reliability  in  a distributed  netv/ork. 

The  decrease  in  the  slope  of  the  network  survivability  curve 
for  sufficiently  small  mobile  networks  is  due  to  the  fact 
that  not  all  nodes  are  of  equal  degree  as  in  distributed  net- 
works. As  a result,  the  loss  of  any  given  node  v/ill  not  have 
as  high  a probability  of  disrupting  the  network. 

It  can  also  be  seen  by  comparing  the  mobi]e  network 
survivability  curves  of  Figures  IV  and  V that  for  corres- 
ponding values  of  the  respective  threshold  reliabilities  in 
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the  eight  node  network,  the  nodes  and  links  affect  the  net- 
work survivability  in  much  the  same  v/ay.  This  leads  the 
author  to  hypothesize  that  if  the  main  consideration  in  a 
sufficiently  small  mobile  network  is  the  cost  of  the  elements, 
then  the  survivability  of  the  network  will  not  be  greatly 
affected  by  choosing  the  lov/cr  cost  element.  For  exairiple , 

suppose  that  the  cost  of  a node  in  the  netv/ork  is  much  greater 

» 

than  the  cost  of  increasing  the  range  of  the  links,  where 
changing  the  restriction  on  the  link  range  is  the  only  feas- 
ible way  to  add  links.  Then  the  hypothesis  is  that  it  would 
be  possible  to  reduce  the  number  of  nodes  and  find  an  in- 
creased range  for  the  links  that  would  realize  a reduction 
in  the  netv/ork  cost  v/ithout  adversely  affecting  the  network 
survivability. 

6.^;  Areas  for  future  development 

Two  aspects  of  mobile  netv/orks  are  viev/ed  as  areas  in 
which  further  v/ork  might  be  useful.  The  first  of  these  would 
involve  the  development  of  the  hypothesis  above  to  determine 
the  effect  on  the  netv/ork  survivability  induced  by  changing 
the  number  of  nodes  and  the  range  of  the  links.  These  tests 
would  also  help  determine  the  maixiraum  size  of  the  mobile  net- 
v/ork for  v/hich  the  survivability  of  the  network  would  not  be 
reduced  when  reducing  the  netv/ork  cost. 

The  second  area  of  possible  study  involves  the  movement 
of  the  nodes  in  the  network  translation  algorithm.  The  algo- 
rithm positions  the  nodes  such  that  the  degree  of  each  node 
is  at  least  as  great  as  the  network  connectivity.  A differer.t 
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Appendix  A 
Program  Listing 


Contained  in  this  appendix  is  a listing  for  the  program 
written  by  this  author  to  implement  the  net^vork  translation 
algorithm.  The  program  v.’as  written  in  FORTRAN  IV  v/ith 
VJATFOR/./ATFI  Y . 
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DIIffiNSIOIM  PST1:(50.4)  ,FSTP(50,4)  ,DISP(50)  ,LINK(2,200)  . 

2 DIST(50 ,50) ,PSLI(200) ,PGK0(50) ,ND£G(50) ,£RAN(50) , 

3 LRAN(200)  ,RANK(50)  ,RANL(200)  .PTEI-U  50 ,4)  ,NI.:0V(50)  , 

4 AR(50,2) ,MVAR(50) ,ISTP(50) ,NMB(50) ,R0K(50) ,NAD(50) 

REAL  NFLG(50)/50*’  ’/. ALTAR/'  '/.MD 

LOGICAL  RANDOM , FLAG , FINAL , FIRST 

DATA  NS,LS/5G,200/ 

DATA  MD , MI  NP , I.’CASE , PC  , RN  , SEP/4 .0,2,20,50.0,4.0,0.5/ 

READ  (5,402)  PCLN,PCND,NRAT 

READ  (5,403)  N,ND,((PSTR(I,J)  ,J=1,ND)  ,I=1,N)., 

2 ((PSTP(I,J) ,J=1.ND) .I=1,N) 

ICASE  = 0 
C 

C CALCULATE  THE  TOTAL  DISPLACEMENT  OF  EACH  NODE. 

C 

DO  10  I = l.N 
NMOV(I)  = 0 

10  DISP(I)  = DISTAN  (NS,I,I,ND,PSTR,PSTP) 

NT  = INT(N*PC*1 .E-2) 

WRITE  (6,501) 

WRITE  (6,502)  (I  ,(PSTR(I  ,J)  ,J^-1  ,ND)  ,(PSTP(I  ,J)  , 

2 J=1 ,ND) ,1=1 ,N) 

WRITE  (6,515)  PCND,PCLN.RN,N,NT,MD,NRAT 
FCTL  = (100.0  - PCLN)*1 .E-2/RN 

FCTN  = (100.0  - PCND)*l.E-2/(N-l .) 

C 

C INSURE  THAT  THE  TERJWNAL  AND  INITIAL  POSITIONS  ARE 

C CONInSCTED. 

C 

CALL  ADJA  (NS,LS,N,ND,RN,PSTP,DIST,LINK,NDEG,r.l) 

T.1NDG  = NDEG(l) 

DO  20  I = 2, N 

20  IF  (NDEG(I).LT.MNDG)  MNDG  = NDEG(I) 

RANDOM  = .FALSE. 

CALL  COMPON  (NS ,N ,M , RANDOM ,PSNO ,PSLI ,RANN ,RANL , 

2 I.I  NK  , IC  OM , N VAR , N V ) 

IF  ( LOOM. ME. 1)  GO  TO  150 

CALL  ADJA  (NS,LS,N,ND,RN,PSTR,DIST,LINK,NDEG,M) 

DO  30  I = 1,N 

30  IF  (NDEG(I) .LT.MNDG)  MNDG  = NDEG(I) 

WRITE  (6,517)  MNDG 

CALL  COMPON  (NS ,N ,M .RANDOM ,PSNO ,PSLI ,RANN ,RANL, 

2 LINK,LCOM,NVAR,NV) 

IF  (LCOM.NE.l)  GO  TO  l60 
FINAL  = .FALSE. 

C 

C BEGIN  THE  MAIN  CALCULATIONS.  CONDUCT  A RANDOM 
C ATTACK  ON  THE  CURRENT  NETWORK  POSITION. 

C 

MO  = NT 

45  WRITE  (6,503) 

WRITE  (6,504)  (l.NFLGd)  ,(PSTR(I,J)  ,J=1,ND)  ,NDEG(I)  , 

2 DISP(I) ,1=1 ,N) 

WRITE  (6,516) 
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36 

DO  47  I ==  1 ,N 

37 

47 

NFLG(I)  = BLAI-^K 

38 

RANDOM  = .TRUE. 

39 

DVD  = 0. 

40 

TLR  = 0. 

4l 

TNR  = 0. 

42 

DO  48  K =.1,M 

43 

I = LINKd.'K) 

4^j 

J = LINK(2,k) 

45 

PSLI(K)  = 1.  - FCTL*DIST(I ,J) 

46 

48 

TLR  =■■  TLR  + PSLI(K) 

4? 

AVLI  = TLR/n 

48 

DO  49  I = 1,N 

49 

PSNO(I)  ==  1.  - FCTN*NDEG(I) 

50 

49 

TNR  = TNR  PSNO(I) 

51 

AVr^O  = TNR/N 

52 

WRITE  (6,521)  AVNO.AVLL 

53 

DO  70  ITER  l.NRAT 

54 

READ  (5,401)  (NRAN(I) ,I=1,N) 

55 

READ  (5,401)  (LRAIULN)  ,LN=l,I.i) 

56 

DO  50  LN  = 1,M 

57 

50 

RANL(LN)  = LRAN(LN)*l.E-2 

58 

DO  60  I = 1,N 

59 

6o 

RANN(I)  = RRAI'I(I)*l.E-2 

6o 

CALL  COWPON  ( NS , K ,M , RANDOM , PSNO , PSLI , RANN , RANL , 

2 LINK,LC0I.1,MVAR,NV) 

61 

70 

IF  (LCCM.SQ.l)  DVD  = DVD  + 1. 

62 

REL  = DVD/NRAT 

63 

WRITE  (6,507)  REL 

64 

IF  (FINAL)  GO  TO  140 

65 

WRITE  (6,508)  I CASE 

66 

I CASE  = I CASE  + 1 

67 

IF  (ICASE.GT.NCASE)  GO  TO  I70 

68 

DO  69  I = 1,N 

69 

69 

ISTP(I)  = 0 

70 

n 

RANDOM  = .FALSE. 

\j 

c 

MOVE  AS  MANY  NODES  AS  POSSIBLE  AND  INSURE  THAT  THE 

c 

p 

RESULTING  NETWORK  IS  CONNECTED. 

71 

DO  71  I = 1,N 

72 

AR(I,1)  = I 

73 

AR(I,2)  = DISP(I) 

74 

DO  71  J = 1,ND 

75 

71 

PTEM(I,J)  = PSTR(I,J) 

76 

CALL  SORT  (NS,H,AR) 

77 

FIRST  = .TRUE. 

78 

DO  77  J = l.MO 

79 

Nr.IOV(J)  = AR(J,1) 

80 

I = NMOV(J) 

81 

77 

CALL  TEMIN  (NS ,N ,ND,MD,I ,DISP( I ) ,ISTP( I ) ,SEP ,PSTR , 

2 PSTP.PTEM) 

82 

80 

CALL  ADJA  (NS,LS,N,ND,RN,PTEM,DIST,LINK,NDEG,M) 
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83 

CALL  COUPON  ( NS , N . M , RANDOM , PSNO , PSLI , RANN , RANL , 

2 LINK, LOOM, N VAR, NV) 

84 

IF  (LCOM.EQ.l)  GO  TO  90 

85 

K = 0 

86 

L = 0 

87 

DO  184  I = l.MO 

88 

DO  181  J - '1,NV 

89 

181 

IF  (NKOY(I) .SQ.NVAR(J))  GO  TO  182 

90 

GO  TO  183 

91 

182 

K = K + 1 

92 

KI^B(K)  ^ NMOV(I) 

93 

GQ  TO  184 

94 

183 

L L + 1 

95 

NOK(L)  = ro.iov(i) 

96 

184 

CONTINUE 

97 

lA  = 0 

98 

IF  (K.EQ.O)  GO  TO  193 

99 

DO  186  I = 1,NV 

100 

DO  185  J = 1,K 

101 

185 

IF  (NVAR(I) .EQ.NM3(J))  GO  TO  186 

102 

lA  = lA  + 1 

103 

NAD(IA)  ==  NVAR(I) 

104 

186 

CONTINUE 

105 

DO  189  I = 1,K 

106 

J = NI.1B(I) 

107 

ISTP(J)  = ISTP(J)  + 1 

108 

IF  {ISTP(J) .GT.MTNP)  ISTP(J)  = MIHP 

109 

189 

CALL  TEMIN  (NS ,N,ND,MD, J,DISP( J) ,ISTP( J) ,SEP,PSTR , 
2 PSTP,PTEI.l) 

110 

GO  TO  195 

111 

193 

DO  194  I = 1,NV 

NAD(I)  = NVAR(I) 

112 

194 

113 

lA  = NV 

114 

195 

IF  (lA.EQ.O  .OR,  L.EQ.O)  GO  TO  80 

115 

DO  187  I = 1,L 

116 

J = NOK(I) 

117 

AR(J,1)  = J 

118 

I87 

AR(J,2)  = NDEG(J) 

119 

CALL  SORT  (NS.N.AR) 

120 

IF  (lA.GT.L)  lA  = L 

121 

DO  190  I = 1,IA 

122 

DO  191  J = l.MO 

123 

191 

IF  (AR(I,1).EQ.NM0V(J))  GO  TO  192 

124 

GO  TO  190 

125 

192 

DO  188  II  = 1,ND 

126 

188 

PTEM(J,II)  = PSTR(J,II) 

K = NAD(I) 

127 

128 

NMOV(J)  = K 

129 

ISTP(K)  = '0 

130 

CALL  TEMIN  (NS,N,ND,MD,K,DISP(K) ,IST?(K) .SEP.PSTR, 
2 PSTP,PTEM) 

131 

190 

CONTINUE 

132 

GO  TO  80 

133 

85 

CALL  ADJA  (NS,LS,N,ND,RN,PTEI.1,DIST,LIN-K,NDEG,M) 
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134 

90 

IF  (FIRST)  IC  = 0 

135 

FIRST  ^ .FALSE. 

136 

DO  93  I = l.N 

IF  (NDEG(I) .LT.MNDG)  GO  TO  92 

137 

93 

138 

GO  TO  99 

139 

92 

DO  94  J = l.KO 

l4o 

94 

IF  (MMOV(  J)  .EQ.I)  GO  TO  96 

141 

II  KO  - IC 

142 

IC  = IC  + 1 

143 

IF  (lI.EQ.O)  GO  TO  99 

144 

NH  KKOV(II) 

145 

NKOV(II)  = I 

146 

DO  95  J = 1 ,ND 

14? 

95 

PTE1.1(NH,J)  = PSTR(NH,J) 

148 

CALL  TEr.IN  (HS,N,KD,MD,I,DISP(I)  ,ISTP(I)  .SEP.PSTR, 
2 PSTP.PTEM) 

149 

GO  TO  85 

150 

96 

ISTP(I)  = ISTP(I)  + 1 

151 

IF  (ISTP(I)  .GT.INT(2.mD-l))  GO  TO  132 

152 

I ONE  = 1 

153 

ONE  = 1. 

154 

CALL  TEI.HN  (NS,N,ND,ONE,I,DISP(l)  .lOIffi.SEP.PTEM, 

2 PSTR.PTEW) 

155 

GO  TO  85  ' 

156 

99 

DO  110  I = l.MO 

157 

J = NI.'IOV(I) 

153 

IF  (J.EQ.O)  GO  TO  120 

159 

DO  100  K - 1,ND 

160 

100 

PSTR(J,K)  = PTEM(J,K) 

161 

DISP(J)  = DISTAN  (NS,J,J,ND,PSTR,PSTP) 

162 

NFLG(J)  ^ ASTAR 

163 

110 

Nr.!OV(l)  = 0 

164 

120 

KO  = 0 

165 

DO  130  I = 1,N 

166 

130 

IF  (DISF(l).GE.l.E-4)  HO  = MO  + 1 

167 

IF  (MO.GE.NT)  MO  = NT 

168 

IF  (MO.EQ.O)  FINAL  = .TRUE. 

169 

GO  TO  45 

170 

132 

J =■-  1 

171 

CALL  TEMIN  (NS ,N ,ND,MD,I ,DISP(I ) , J ,SEP ,PSTR ,PSTP , 

2 PTEH) 

172 

DO  135  K = l.N 

173 

IF  (I.EQ.K)  GO  TO  135 

174 

AR(K,1)  = K 

175 

AR(K,2)  = DISTAN  (NS,I ,K,ND,PTEM,PTEM) 

176 

135 

CONTINUE 

177 

CALL  SORT  (NS,N,AR) 

178 

Ml  = N - MNDG 

179 

K = AR(M1,1) 

180 

DOL  = AR(M1,2)  - MD  + O.O5 

181 

DO  138  L = 1,ND 

182 

138 

PSTR(I,L)  = PTEM(K,L) 

183 

J = 0 
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184  CALL  mil N (NS,N,ND,D0L,I,AR(l.ll,2j  ,J,SEP,PTEM, 

2  PSTR.PTEM) 

185  GO  TO  85 

186  140  WRITE  (6,510)  ICASE 

18?  WRITE  (6,601) 

188  GO  TO  200 

189  150  WRITE  (6,511) 

190  GO  TO  200 

191  160  WRITE  (6,512) 

192  GO  ‘ixi  200 

193  170  WRITE  (6,513) 

C 

Q PORiviatS****************** 

194  401  FORMAT  (4012) 

195  402  FORMAT  ( 2( 6X , F? . 2 ) , lOX , I4) 

196  403  FORMAT  (2I7/(10X,2F15.6)  ) 

C 

Q p0pj,'iATS^'********‘-^*^'”'‘*'*‘‘'*'*'‘* 

197  501  FOPJaT  ( *1 '//1 5X,  ' INITIAL  AND  FINAL  POSITIONS  OF', 

2 ■•STATICNS'//24X,  'INITIAL*  ,11X,  'FINAL* ) 

198  502  FORMAT  (14X-,I5,2F7.1  ,3X,2F7.1) 

199  503  FORMAT  ( '1'//15X, 'NODE  X Y DEG  DISTANCE'//) 

200  504  FORMAT  ( l6X ,12 , A2 , 2F5 . 1 ,l6 , Fll . 3) 

201  507  FORMAT  ( ///l OX , 'PROBABILITY  TFJ^T  NETWORK  SURVIVES  ', 

2 'RANDOM  ATTACK: • ,F6. 3)  . 

202  508  FORMAT  (/1 5X, 'TIME  INTERVAL* ,14) 

203  510  FORMAT  (//l OX, 'TERMINAL  POSITIONS  REACHED  IN  TIME  '. 

2 'INTERVAL' ,15) 

204  511  FORMAT  (//I5X ,' DESIRED  TERMINAL  POSITIONS  ARF  NOT  ', 

2 'CONNECTED  USING  GIVEN  RANGE  RESTRICTION') 

205  512  FORMAT  ( //I 5X, 'DESIRED  INITIAL  POSITIONS  ARE  NOT  ', 

2 'CONNECTED  USING  GIVEN  RANGE  RESTRICTION') 

206  513  FORMAT  (//I OX, 'MAXIMUM  NUMBER  OF  Tllffi  INTERVALS  ', 

2 'EXCEEDED  'WITHOUT  REACHING  DESIRED  TERMINAL  POSI ' , 

3 'TIONS') 

207  515  FORMAT  (///18X, 'MINIMUM  NODE  RELIABILITY: ' ,F6.1 , 

• 2 ’ ' /1 8X, 'MINIMUM  LINK  RELIABILITY:  ' ,F6.1 , ' ^'/24X, 

3 'MAXIMUI.l  LINK  RANGE:  ' ,F7.2/21X, 'TOTAL  NUI.’.BER  OF  ', 

4 'NODES: ' ,I4/20X, 'NUMBER  OF  NODES  MOVING '/25X , 'PER  ', 

5 'TIME  INTERVAL: 'I4/17X, 'MAXIMUIii  NODE  DISPLACEMENT'/ 

6 25X,'PER  TIME  INTERVAL: ' ,F7.2/18X, 'NUMBER  OF  RAN', 

7 'DOM  ATTACKS '/25X, 'PER  TIME  INTERVAL: ' ,I4) 

208  516  FORMAT  (//19X,'*  NODE  MO\'ED  DURING  CURRENT  TIME  ', 

2 'INTERVAL') 

209  517  FORMAT  (/I 3X, 'MINIMUM  DEGREE  ( CONNECTIVITY) */20X , 

2 'ALLOWED  ON  MOVED  NODES: ',I4) 

210  521  FORMAT  (//1 5X, 'AVERAGE  NODE  RELIABILITY: ' ,F8.3/15X, 

2 'AVERAGE  LINK  RELIABILITY :', F8 . 3) 

C 

211  200  STOP 

212  END 
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213  FUNCTIOri  DISTAN  (NSIZE.I , J,ND,PBEG  ,P2ND) 

C 

C THIS  ROUTiriE  CALCULATES  AND  RETURNS  THE  DISTANCE 
C BSTVSEN  POINT  PBSG(I,K)  AND  PEND(J,K)  FOR 

C K=1,...,HD. 

C 


214 

DIMENSION 

PBEG ( NSIZE , ND) , PEND ( MSI ZE , MD) 

215 

DIS  = 0. 

216 

DO  10  K = 

1,ND 

217 

10  DIS  = DIS 

+ (PEr!D(J,K)  - PBEG(I,K))**2 

218 

DISTAN  = : 

:SQRT(DIS) 

219 

RETURN 

220 

END 

k 


221 


SUBROUTI I;E  ADJA  ( NS  , r.r.  N'  , ND  , RN  , POSI  , di st  , link  , i deg  , ?0 

c 

C THIS  ROUTINE  DETERMINES  VJHICII  NODES  ARE  ADJACENT. 

C ' 

G RETURNED  VALUES: 

C . LINK(2,]0:  TERMINAL  NODES  OF  THE  LI InKS. 

C IDEG(N)  : DEGREES  OF  THE  NODES. 

C K : ACTUAI.  NUMBER  OF  LINKS  IN  NETVJORK. 

C 


222 

DIMENSION 

DIST(NS,N) ,P03I(NS,ND) ,LINK(2,M) ,IDEG(N) 

223 

K = 0 

2Zk 

DO  10  I = 

l.N 

225 

10  IDEG(I)  = 

0 

226 

1 

li 

1 

227 

DO  20  I = 

1 ,NI/il 

228 

IPl  = I + 

1 

229 

DO  20  J = 

IPl.N 

230 

DISTd,  J) 

= DISTAN  (NS.I.J.ND.POSI.POSI) 

231 

IF  (DISTd  ,J)  .GT.RN)  GO  TO  20 

232 

K = K + 1 

233 

LII^;K(1  ,K) 

= I 

234 

LINK(2,K) 

= J 

235 

IDSGd)  = 

IDEGd)  + 1 

236 

IDEG(J)  = 

IDEG(J)  +1 

237 

20  CONTINUE 

238 

RETURN 

239 

END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTI  im  COI.1POM  ( NM  , N , M , RANDOM  , PSN , PSL , RAN  , RAL , 

2 LlNK,LCOI.i,NV,K) 

THIS  ROUTINE  DETERMINES  THE  MUI'.IBER  OF  COMPONENTS  IN 
THE  NETWORK  BOTH  BEFORE  AND  AFTER  A RANDOM  ATTACK. 
WHEN  RANDOM  =-■  .FALSE.  THE  ROUTINE  CHECKS  THE  CON- 
NECTIVITY OF  A DETERMINISTIC  NETl.DRK.  ’.HiEN  RANDOM 
.TRUE.  THE  CONNECTIVITY  IS  CHECKED  FOR  A PROBABI- 
LISTIC NETWORK  WHOSE  NODES  AND  LINKS  FAIL  RANDOMLY. 

RETURIffiD  VALUES : 

■ LOOM  : NUMBER  OF  COMPONENTS  IN  THE  NETWORK. 
NV(K):  MODES  NOT  IN  ThB:  MAIN  COMPONENT. 


2^1-1 

DIMENSION  NSUR(50) ,NKUM( 50) ,NC0M(50) ,NrffiX(50) . 

2 MI.AS(50)  ,NFIR(50)  ,PSN{N)  ,PSL(M)  ,RAN(N)  ,RAL(M)  , 

3 LINK( 2 ,M ) . AR ( 50 . 2 ) , NV( N ) 

242 

LOGICAL  RANDOM 

243 

DO  10  I = l.N 

244 

NCOM(I)  = I 

245 

10 

NNEX(I)  = 0. 

246 

IF  (RANDOM)  GO  TO  30 

24? 

DO  20  I = l.N 

248 

NNTH.ld)  = 1 

249 

NFIR(I)  = I 

250 

20 

NLAS(I)  = I 

251 

GO  TO  55 

9 <9 

30 

DO  50  I = 1,N 

253 

IF  (RAN(I) .LT.PSN(I) ) GO  TO  4o 

254 

NSUR(I)  = 0 * 

255 

NKUI.I(I)  = 0 

256 

NFIR(I)  = 0 

257 

NLAS(I)  = 0 . • 

258 

, GO  TO  50 

259 

40 

NSUR(I)  = 1 

260 

NNUM(I)  =1 

261 

NFIR(I)  = I 

262 

NLAS(I)  = I 

263 

50 

CONTI  runs 

264 

55 

DO  100  LN  = 1,M 

265 

10  = LINKd.LN) 

266 

JO  = LINK(2,LN) 

267 

IF  ( .NOT. RANDOM)  GO  TO  60 

263 

IF  (RAL(LN) .GE.PSL(IN))  GO  TO  100 

269 

IF  (NSURdO)  .EQ.O  .OR.  NSUR(  JO)  .EQ.  O)  GO  TO  100 

270 

60 

IF  (NCOMdO)  .EQ.NCOM(  JO) ) GO  TO  100 

271 

IF  (HNUr/IdO)  .GE.roiUT>'l(  JO) ) GO  TO  70 

272 

KO  = JO 

273 

JO  = 10 

274 

10  = KO 

275 

70 

I = NCOMdO) 

276 

J = NCOM(JO) 

277 

IJ  = NLASd) 

278 

JJ  = NPIR(J) 
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279 

NNEX(IJ)  - JJ 

280 

KLAS(I)  = KLAS(J) 

281 

KFIR(J)  = 0 

282 

NMS(J)  = 0 

283 

NNUMdO)  = NNUr.l(IO)  + NNUI.1(, 

284 

II  = NFIR(I) 

285 

80 

NNUlKll)  = HNUIl(IO) 

286 

II  = MliEXdl) 

287 

IF  (II.NE.O)  GO  TO  80 

288 

90 

NGOr.l(JJ)  - I 

289 

JJ  = NMEX(JJ) 

290 

IF  (JJ.NE.O)  GO  TO  90 

291 

100 

CONTINUE 

292 

LCOM  = 0 

293 

DO  110  1 = 1 ,N 

294 

IF  (NFIRd)  .EQ.O)  GO  TO  110 

295 

LCOM  = LOOM  + 1 

296 

AR(LC0M,1)  = I 

297 

AR(LC0M,2)  = NNUMd). 

298 

110 

CONTINUE 

299 

IF  (RANDOM)  -RETURN 

300 

IF  (LCOM.EQ.l)  RETURN 

301 

CAli  SORT  (MM.LCOm.AR) 

302 

J = 0 

303 

K = 1 

304 

DO  130  I = 2, LCOM 

305 

J = J + K 

3u6 

K = ARd,2) 

307 

K = J + K - 1 

308 

IJ  = ARd.l) 

309 

IJ  = NFIRd  J) 

310 

DO  130  II  = J,K 

311 

NVdl)  = IJ 

312 

130 

IJ  = NNEXCIJ)  i 

313 

RETURN 

314 

END 
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315 

SU5R0UT1  M-i  ^ 

C 

C 

THIS  RO;  l i. 

c 

DESCEHDlHl.  i-i.. 

c 

c 

RETURHi;D  VA, 

c 

ARC.  .. 

c 

c 

c 

316 

DiMEiiSici':  a:u 

317 

LOGICAL  j l/.. 

318 

LIU  N - 1 

319 

20 

FLAG  = 

320 

DO  60  I 1 . 

321 

J ^ I 1 

322 

IF  (AR(I.2). 

323 

FLAG  = . T2’  : . 

324 

DO  40  K ^ ' 

325 

HOLD  = a;;(: 

326 

AR(I,K)  Ak(. 

327 

40 

AR(J.K)  HCLL 

328 

60 

COHTII.'UE 

329 

IF  (FLAG)  Cl  . 

330 

RETURN 

331 

END 
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332  SUBROUTIKIi  TLIvlIK'  ( rU'.l , N , HD , UK  , I , US  , IG  , SEP  , PS  , PF , TP ) 

c 

C THIS  ROUTiriE  CALCULATES  THE  TEMPORARY  POSITION  FOR 

C THE  ITH  NODE  AND  INSURES  THAT  IT  IS  MO  CLOSER  THAN 

C SEP  TO  ANY  OTHER  NODE. 

C 

C RETURNED  VALUES : 

C TP(I,NU):  TEMPORARY  POSITION  OF  THE  ITH  NODE. 

C 


333 

DIMENSION  PS(NM.ND)  ,PF(NM,ND)  ,TP(NT1.ND)  ,MT(50) 

334 

IF  (DS.LE.DM)  DS  = DM 

335 

DO  100  J = 1,ND 

336 

100 

TP(I,J)  =-■  PS(I,J)  + {PF(I,J)-PS(I,J))-»DI.1/(DS*2^-*IG) 

337 

DO  5 K = 1,N 

338 

5 

MT(K)  = 0 

339 

10 

DO  90  J l.N  ■■ 

340 

IF  (I.SQ.J)  GO  TO  90 

341 

D = DISTAN  (NM.I  ,.T,ND,TP,TP) 

342 

IF  (D.GE.SEP)  GO  TO  90 

343 

IF  (D.LT.l.E-5)  GO  TO  30 

344 

DO  20  K = 1,ND 

345 

20 

TP(I.K)  = TP(I.K)-(TP(J,K)-TP(I,K))*(8EP-D)/D 

346 

GO  TO  50 

347 

30 

DO  40  K = 1,ND 

TP(I,K)  = TP(I,K)-(TP(J,K)-TP(I,K))^-SEP 

348 

40 

349 

50 

D = DISTAN  (1M,I,I,ND,PS,TP) 

350 

IF  (D.LE.DM)  GO  TO  90 

351 

IF  (MT(J).Sw.6)  go  to  70 

352 

MT(J)  = MT(J)  + 1 

353 

DO  60  K = 1,ND 

354 

60 

TP(I,K)  = TP(I,K)-(PS(I,K)-TP(I,K))*(DM-D)/D 

355 

70 

IF  (MT( J) .EQ.IO)  RETURN 

356 

KT(J)  = I.IT(J)  + 1 

357 

DO  80  K = 1,ND 

358 

80 

TP(I,K)  = TP(I,K)-(P3(I,K)~TP(I,K))*DM-D-SEP)/D 

359 

. GO  TO  10 

360 

90 

continue: 

361 

RETURN 

362 

END 
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